Astronomy & Astrophysics manuscript no. bgrovesIR 


February 5, 2008 


(DOI: will be inserted by hand later) 





The Infrared Emission from the 
Narrow Line Region 

Brent Groves 1 , Michael Dopita 2 , & Ralph Sutherland 2 



O 
O 

(N 

< 



> 

in 
o 

oo 

o 

\o 
o 

6 



X 



1 Max Planck Institute for Astrophysics, 
Karl-Schwarzschild str. 1, Garching, 85741 Germany 
e-mail: brent@mpa-garching.mpg . de 

2 Research School of Astronomy & Astrophysics, The Australian National University, 
Cotter Road, Weston Creek, ACT 261 1, Australia 

Received jdate^, / Accepted [date^ 

Abstract. We present models for the mid- and far- infrared emission from the Narrow Line Region (NLR) of Active Galactic 
Nuclei (AGN). Using the MAPPINGS in code we explore the effect of typical NLR parameters on the spectral characteristics of 
the IR emission. These include useful IR emission line ratio diagnostic diagrams for the determination of these parameters, as 
well as Star formation-AGN mixing diagnostics. We also examine emission line to continuum correlations which would assist 
in separating the IR emission arising from the NLR from that coming from the inner torus. We find for AGN like NGC 1068 
and NGC 4151 that the NLR only contributes ~10% to the total IRAS 25 /jm flux, and that other components such as a dusty 
torus are necessary to explain the total AGN IR emission. 



1. Introduction 

Active galactic nuclei (AGN) are well known to emit strongly 
at all wavelengths, from gamma-rays and X-rays to the in- 
frared and radio. Their emission characteristics at each wave- 
length gives us insight into the central engine of these galaxies, 
whether observed directly, such as at X-ray wavelengths, or in- 
directly through reprocessed emission, such as by study of the 
emission lines from the narrow line region (NLR). 

One area that has seen recent rapid advances is studies of 
the mid- and far-infrared characteristics of AGN, with space 
telescopes such as IRAS, ISO and Spitzer making this wave- 
length regime readily accessible. As the continuum in this 
wavelength region is dominated by warm and hot dust re- 
emission, it is particularly useful in testing the paradigm of the 
unified model, which argues for a dense equatorial band of dust 
surro unding the black hole central engine ( Antonucci & Miller, 
1985). 



Various attempts have been made to model the mid- and 
far-IR emission. These range from simple empi rical fits of dust 
grey- b ody emission to t he observations fedelson & Malkarl 
1986; iBlain et all 120031) . to more complicated distributions 
of dust which include the stoc hastic heating of dust by 
an A GN power-law continuum JPullemond & van Bemmeli 
2005). With the prevalence of the unified model, most anal- 
yses of the IR from AGN have concentrated on the emis- 



source (Pier & Krolik, 1992, 1993; Granato & Danese, 


Efstathiou & Rowan-Robinson, 


1995 




Granato et al., 



Il997llvan Bemmel & Dullemondtl2003HNenkova et all E002). 
However, there have bee n some authors who have co nsidered 
different geometries (i.e. ISiebenmorgen et all 120041) or have 
considered i n more detail the co llisional heating and destruc- 
tion of dust (Contin i et allEooH) . 

Often ignored in the considerations of AGN IR emission is 
the contribution from the narrow line region (NLR). Th is re- 
gion is known to contain dust (e.g. Radoms ki et all l2003) and, 
like the dusty torus, will be h eated by the cent ral engine and re- 
emit at IR wavelengths jNetzer & LaoJ.fl993h . With lower den- 
sities than found in the central torus, the NLR also emits strong 
forbidden lines such as the [O ra] 500.7 nm line. Such lines are 
often one of the main diagnostics for the presence of an AGN, 
and also provide important constraints on the form and inten- 
sity of the radiation field which excites them. They also provide 
insight into the physical conditions such as density or pressure, 
and constraints on the metallicity of the circum-nuclear ISM. 

As AGN are often heavily dust-embedded, much of the line 
emission from the near nuclear parts of the NLR may be to- 
tally obscured at optical wavelengths. To penetrate this circum- 
nuclear dust while retaining adequate spatial resolution to in- 
vestigate the structure of the dense parts of the NLR, we need 
to observe at mid-IR wavelengths. Instruments such as the IRS 
instrument of the Spitzer Space Telescope are finally provid- 
ing the data needed to study this region. In order to interpret 
this data we need both line and continuum diagnostics in the 
mid-IR. 

Here, we will explore the expected characteristics of the in- 
frared emission from the NLR, examine the dependence of this 



2 



B. Groves, M. Dopita, & R. Sutherland: The IR Emission from the NLR 



emission on the parameters that define the NLR, and investigate 
possible diagnostics that can be used to define these parameters 
and so give us insight into the NLR and the AGN as a whole. 

2. Modelling the Dust in Narrow Line Regions 

The effects of dust on the structure of the NLR and its influ- 
ence of the emission lin e spectrum has be en discussed in detail 
in previous papers ddroves et alll2004albh . Here we extend this 
work using a more realistic distribution of dust and considering 
in addition the stochastic heating of the dust to investigate the 
nature of the far-IR dust re-emission spectrum and its relation- 
ship to the line spectrum of the NLR. 

The existence of dust around AGN is well accepted. Indeed, 
it forms an inte gral component o f the unification scenario of 
active galaxies dAntonuccH Il992i) . That dust specifically ex- 
ists wit hin the NLR of AGN is also suggested by theoretical 
models dNetzer & Laorl Il993l bopita et all 120021) and is well 
suppor ted by recent high r esolution IR observations of nearby 
AGN dAlloin et all bOOOi iBock et all l2000HPackham et all 
l2005UGalliano et allbOOfllMason et alll2006l) . 

Both the composition and size distribution of the dust in the 

NLR are still uncertain. Current observations (Gaskell et al., 

■ 1 

2004) suggest that the properties of dust in the vicinity of AGN 

may differ to that found in our Galaxy. This would be unsur- 
prising considering the gas pressure, UV radiation field and 
dynamics of this environment are all much more extreme than 
those found in the solar vicinity. Polycyclic aromatic hydrocar- 
bons (PAHs), while common in our Galaxy and are observed 
in starforming galaxies, are not expected in the excited regions 
of active galaxies. This is because of the strong EUV and X- 
ray emission from the AGN, which are thought to effecti vely 
destroy PAHs though photo-dissociation dVoillll99 llll992h . In 
the cases where PAH emission is seen in the nuclear regions of 
AGN, the emission is thought to arise from nuclear starbursts. 
The starbursts supply the Far-UV needed to excite the PAH 
emission, but would be shielded from the X-ray e mission of 
the AGN by intervening dusty , obscuring material dlmanishil 
l2003Ulmanishi & Wadalbool . 

The lack of these planar linked benzene-ring molecules 
may also extend to small grains, with the smallest graphite and 
silicate grains possibly being simila rly destroyed by the st rong 
radiation field. As first suggested bv lLaor & DraineN 19931) . the 
idea that larger grains dominate the dust in AGN relative to 
that found in the Galaxy has been supported by observations 
such as the recent work of Gaskel l et al.l (I2004I) . which demon- 
strated the muc h flatter attenuation curves in AGN in the UV. 
iMai olino etalJfeOOIalhl) explored the observational evidence 
for different dust properties within AGN and came to the con- 
clusion that the most favorable scenario was one where large 
dust grains dominated the distribution of AGN dust. However, 
flatter attenu ation curves may also arise in a h ighly clumpy dust 
distribution dFischera & DopitaH2004[ l 2005l) . 

The composition of dust within AGN is also poorly con- 
strained by observation. Metals appear to be depleted within 
the NLR gas, although the magnitude of this depletion is un- 
certain. In terms of extinction, the standard silicate + graphite 
dust appears to reproduce the observations reasonably well 



(Laor & DraineL 1 19931 IMaiolino et all l2001albl) . although ad- 
ditional dust component s, such as diamond dust have also been 
suggested dBinette et alll2005i) . 

Using this information we have tried to construct a sim- 
ple yet physically plausible model for the dust in the NLR 
clouds. We assume two types of dust: graphitic dust and 
silicaceo us dust. The o p tical d a ta for these are taken from 
work of iDraine & Led dl984l) . Laor & Drainel dl993h and 
IWeingartner & Draineld200lh . 

The grain size distribution of both types of dust is taken to 
be a power-law, with the index arising from grain shattering. 
Grain shattering has been shown to lead naturally to the for- 
mation of a power- law size distribution of gr ains with index 
a ~ 3.3 djones. Tielens & H ollenbachl 1 19961 ). This is within 
the bounds of the standard Mathis, Ru mpl. & Nordsieckld 19771) 
value of a ~ 3.5. Such a grain size distribution can also hold 
only between certain limits in size, the minimum size being 
determined by natural destruction processes such as photode- 
struction, and the upper limit being determined by limits on the 
growth by condensation and sticking. These may occur within 
the AGN or previous to the dusty gas entering the energetic re- 
gion. To capture these elements of the physics, we have adopted 
a modified grain shattering profile with the form; 

e -(a/fl min )- 3 

dN(a)/da = ka a — -r. (1) 

This creates a smooth exponential cut-off in terms of the 
grain mass at both the minimum and maximum grain size of the 
distribution. This smooth cut removes any edge effects in either 
the emission or extinction by dust that arise due to the sharp 
cut-offs in the other distributions. To implement the larger grain 
dominated distribution discussed above, we have assumed that 
the smallest grains have been destroyed due to the energetic 
environment, and used a larger minimum grain size relative to 
Galactic of a^m = 0.01 fim. 

The constant k is determined by the total dust-to-gas 
mass ratio, which is determined by the depletion of the 
heavy elements onto dust. Our models use a solar abun- 
dance set (as modified by the latest a b undance determi - 
nations lAllende Prieto etafl d200ll 120021) : lAsplundl l|2000); 
lAsplund et all fcOOOk The depletion factors are those used 
by iDopita et al.l d200fi) for starburst and active galaxy pho- 
toionization modeling and are similar to those found by 
Ijenkins. Jura & LoewensteirJ dl987l) and ISavage & Sem bach 
dl996h in the local ISM using the UV absorption lines to probe 
various local lines of sight. 

Within the MAPPINGS inr code, the dust grain size distri- 
bution is divided into 80 bins spaced logarithmically between 
0.001-10 fim. The number of grains of each type in each bin 
is then determined by equation ^ The absorption, scattering 
and photoelectric heating is calculated for each size bin. This 
is then used to establish the temperature probability distribution 
for computation of the FIR re-emission spectrum with quantum 
fluctuations. 
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3. Calculating the IR Emission from Dust 

In energetically active regions such as found around starbursts 
and AGN it is essential to take stochastic quantum heating of 
dust grains properly into account. In particular, the temperature 
fluctuations of grains serves to provide a population of small 
grains which are hotter than their equilibrium temperature and 
thus which show an excess in their mid-IR emission. 

Temperature fluctuations and stochastic heating processes 
have been considered by se veral authors since G reenber g first 
proposed these ideas (e.g. fpulevi fl973t iGreenberg & Hongl 
1 19741 IPurcellL 1 19761: IPwekL Il986t) . The most detailed work 
on dust and temperature fluctuat ions has been a se r ies of 
papers by Draine & c o-authors JPraine & Andersoil Il985t 
IGuhatha kurta & Draine[|l9 89 | ) culminating in a set of papers 
with Li JPraine & LiLl200lULi & Drainell200ll) . 

The FIR emission treatment used in the M APPING S hit 
code i s based on the algorithms of IGuhathakurta & Drainel 
Jl989l) (GD89) and lDraine & UfeOOll) (PLOD. The code pro- 
vides an optimized solution of their algorithms to determine 
the dust grain temperature distributions for each grain size ac- 
cording to the both the strength and detailed spectrum of the 
local radiation field. The code then integrates the resultant FIR 
emission of the ensemble of dust grains to provide the local re- 
emission spectrum, which is itself in tegrated throug h the model 
in the outward-only approximation JGroveslEool . 

In general, the average temperature of dust grains decreases 
from the front of a photoionized cloud exposed to the radiation 
to the back as expected. However the average dust tempera- 
ture is not an exact representation of the grains, with the largest 
grains at lower temperatures, then the temperature of the grains 
increasing with decreasing grain size. At a certain size stochas- 
tic heating becomes important and the probability distribution 
of the dust peaks at a temperature below the total average dust 
temperature, but with a tail extending to very high tempera- 
tures. 

4. Narrow Line Region Models 

The physical details of the dusty, radiation-pressure dominated 
photoionization models used here to model the IR continuum 
of the narrow line region clouds have been discussed in depth 
in previous papers ( Gro ves et all l2004albl) . Here we only pro- 
vide a brief overview of the different input parameters and the 
parameter space explored. 

4.1. Input Ionizing Spectrum 

A simple power-law radiation field, with an index around ~ 
- 1 .4 and extending from the FU V to X-ray, is able to reproduce 
the dominant features observed in NLRs (see e.g. Osterbrock, 
Il989i) . However, when the full spectral energy distribution of 
the NLR is being considered a more physical representation 
of the ionizing spectrum is required. The main problem here 
is that the intrinsic spectrum emitted by the NLR is generally 
heavily obscured by the circum-nuclear interstellar matter in 
the active galaxy before it reaches the observer. This is espe- 
cially true in the extreme-UV. 



Our understanding of this region of the spectrum relies pre- 
dominantly on photoionizati on theory and on the theory of ac- 
cretion disk models (see eg. lAlexander et al1ll999t E"000). This 
method gives insight into the spectrum, but is somewhat depen- 
dent upon the modelling used and the other model parameters. 

However, rec ent work from IScott et a"fl J2004I) and 
IShang etalJll2005l) looked at QSOs beyond the Lyman limit us- 
ing space telescopes, exploring not only the peak of the big blue 
bump of the QSO spectra but also the slope beyond this into 
the EUV. These new observations help to constrain the general 
shape of the AGN spectrum and thus the photoionization and 
accretion disk models. 

We have chosen here to base our in put ionizing spec - 
trum predominantly on the o bservations of Elvis et alJ i 19 94). 
while still taking note of fhe lScott etail J2004I) and lShang et alJ 
d2005l) data. To fit the data we us e a combination of two power- 
laws with exponential cut-offs dNagao et all l200ll iFerlandt 
1996), 

fy = v^Kph^W-nT) (2 ) 

+ av°* exp (-^)ex P (-^). 

The parameters are chosen to give the most plausible fit, con- 
centrating on the ionizing spectrum. The final spectrum, shown 
in figure^ is similar to the 3-part broken power-law often used 
in photoionization models. The first parameter, q-euv, is defined 
as the EUV power-law index, and hence controls the dominant 
part of the spectrum. It is chosen here to be ffEuv = -1-75, 
which gives a smooth fit from the peak of the big blue bump 
(BBB) and the soft X-ray. This fit is made smooth by the UV 
cut-off, with kTuy = 120 eV. The parameter £7bbb defines the 
peak of the BBB, and is set to ATr rr =7.0 eV. which gives a 
peak around 1 10.0 nm, as found by Shan g et al 1 il2004 . 

The second component of equation|2]represents the X-ray 
part of the AGN spectrum, which has an index of ax = -0.85 
and an upper cut-off simply chosen to be kTx = 10 5 eV to 
prevent errors. The parameter a sets the scaling between the 
two components of the AGN spectrum, and is se t to a — 0.0055, 
which gives not only a good fit to the E lvis et all i 19941) data but 

also gives an optical-X-ray slope 1 of index ao-x 1«4, fitting 

in with previous photoionization models. 

4.2. Narrow Line Region Metallicity and Dust 
Depletion 

In previous work JGroves et alll2004albl) we found that the ob- 
served line ratios for NLR are consistent with a heavy element 
abundance of between solar and twice solar metallicity. Here, 
for simplicity, we compute mod els wi t h sola r metallicity, us- 
ing the abundances from iDopita et all ( 120051 with references 
therein). For the depletion factors, we follow our previous work 
and use those found in the local ISM due to our ignorance about 
what these factors are within the AGN environment. The solar 
abundances and depletion factors used are given in tabled 

1 Defined as in equation 1 from lNaeao et alJ 1200 ll) 

2 All abundances are logarithmic with respect to Hydrogen 

3 Depletion given as log(X/H) gas - log(Z/H) ISM 
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Fig. 1. The input model ionizing spectrum (scaled) from equa- 
tion|2 

Table 1. Solar abundance & metallicity scaling 



Element Abundance 2 Depletion 3 

H O000 000 

He -0.987 0.00 

C -3.61 -0.30 

N -4.20 -0.22 

O -3.31 -0.22 

Ne -3.92 0.00 

Na -5.68 -0.60 

Mg -4.42 -0.70 

Al -5.51 -1.70 

Si -4.49 -1.00 

S -4.80 0.00 

CI -6.72 -0.30 

Ar -5.60 0.00 

Ca -5.65 -2.52 

Fe -4.54 -2.00 

Ni -5.75 -1.40 



4.3. Choice of Variables 

The density and ionization structure of dusty, radiation- 
pressure dominated photoionization model s used to represen t 
the NLR has been described in detail in lGroves et alJ J2004a). 
Here, we follow this prescription in most respects. However, 
we here define the initial parameters in a slightly different man- 
ner in order to provide the clearest distinction in distinguishing 
the effects of each parameter upon the infrared emission. 

As discussed in the previous work, it is the total pressure 
of the system which defines the density structure of the ion- 
ized plasma. The total pressure is the sum of both the local gas 
pressure and the incident radiation pressure and can be set as 
a parameter through setting both of these quantities. In a static 
cloud, assumed here, the (negative) gradient of the local radia- 
tion pressure must match the (positive) gradient of the gas pres- 
sure. The effective radiation pressure imparted as gas pressure 
at any point in the cloud depends upon both the local ionizing 
spectrum (both shape and flux density) and the optical depth (a 
function of both the column density and the opacity of the gas). 



As our models have the same ionizing spectral shape, the only 
parameter at a given optical depth is the incident flux density 
of the ionizing radiation. 

Thus at a given optical depth, such as near the ioniza- 
tion front (where the dominant ionizing radiation has been ab- 
sorbed), we have: 

^tot — ^gas ^rad 
= ^Ogas + ah)- 

fogas is the initial gas pressure of the NLR cloud, a is a constant 
which includes the effect of spectral shape and gas and dust 
opacity, and Iq is the incident ionizing flux density upon the 
NLR cloud. 

This results in the same density structu re as shown in 
iDopita et all d2002l) and ldroves et alJ J2004ah . however with a 
much better defined input. The definition of density structure is 
important as it allows us to distinguish which variations in the 
IR emission arise from changes in the ionizing flux or the NLR 
density structure. 

4.4. Parameter Space of Models 

Within this work we explore three different total pressures, 
P tot /k =s 10 6 , 10 7 , and 10 8 KcnT 3 , which simulate [Sn] den- 
sities of approximately 10 2 , 10 3 , and 10 4 crrT 3 respectively. 

For each pressure regime we have investigated a range of 
ionizing flux densities. This range has been chosen to give a 
reasonable representation of the typical ionization conditions 
found within the NLR. The values of the input flux densi- 
ties and corresponding initial gas pressures are given in ta- 
ble |2 The total flux density has been normalised by the fac- 
tor 3.8098((P tot //t)/10 6 )ergcm^ 2 s- 1 (or for ionizing flux den- 
sity by 2.0864((P tot //fc)/10 6 ) ergcirrV 1 ). This factor results in 
a ionizing photon density centered around S„/c = 1 for the 
Ptot/k = 10 6 models, giving a range of ionization parameters 
centered approximately around 10~ 2 . 

Traditionally, the ionization parameter t/ 4 has been used 
to describe photoionization equilibrium conditions. However 
in isobaric systems this parameter is harder to define. So we 
present here two parameters for comparison. The standard ion- 
ization parameter as a reference to previous work, and the S 
parameter, 

Pgas nkTc ' 

a more suitable paramete r to describe the io nization equilib- 
rium in isobaric systems (K rolik et allll98ll) . To convert Hi on 
to include the total incident flux density rather than just the ion- 
izing flux, multiply by 1.826. 

As the initial density is uncertain in isobaric models, we 
define the ionization parameter in terms of a normalised initial 
ionization parameter, Uq. This parameter assumes a tempera- 
ture of 10 4 K to obtain the initial density from the input gas 
pressure. As the initial temperature usually exceeds 10 4 K, and 

4 The ionization parameter U is a dimensionless measure of the 
number density of ionizing photons (S ,), over the gas density; U = 
SJn H c. 
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Table 2. Input parameters for NLR models. The total incident 
flux density / is scaled by 3.8098((P lot /k)/l0 6 )eTgcm- 2 s- 1 
(see text), while Pog^/k 1S scaled by Q.l(P tot /k) 



h 




Ho 


log((/ ) 


4.0 


1.00 


20.16 


-0.40 


3.0 


3.20 


4.73 


-1.03 


2.0 


5.40 


1.87 


-1.43 


1.0 


7.60 


0.66 


-1.88 


0.5 


8.70 


0.29 


-2.24 


0.25 


9.25 


0.14 


-2.57 


0.1 


9.58 


0.05 


-2.98 



the initial temperature increases with ionization parameter, the 
normalised initial ionization parameter underestimates the true 
initial ionization parameter. This becomes increasingly so for 
higher values. The value of both Uo and Ho for each model are 
given in table |2 

The models are truncated at a column density of 
log( /V(Hi)) _ 21.5, a reaso nable estimate for the NLR clouds 
fe.p. lCrenshaw et alll2003l) . which means that we are compar- 
ing like absorbing columns when determining the IR emission. 

5. Narrow Line Region SEDs 

The resulting Spectral Energy Distributions (SEDs) for the 
10 7 model are shown in figure |3 The figure shows four SEDs 
corresponding to S = 20.16, 1.87,0.29,0.05. 

In this figure we see several features which arise due to 
the combination of the changing ionizing flux density and the 
changing ionization parameter. The offset between the SEDs 
arises due to the difference in input flux density. A greater in- 
cident ionizing flux density results in a greater nebular suface 
brightness. However, the correlation between the resulting neb- 
ula surface brightness or emission measure and incident ioniz- 
ing flux density is not one to one, due to the varying ratios 
of dust absorbtion compared with absorption in the gaseous 
phase. 

The ionization parameter is one of t he main det erminants 
of a nebula emission spectrum (see e.g. Donita & Sutherland. 
120031 and also re lates to the dust dominance of ionizing 
photon absorption (Dop ita et all [2002). This parameter has a 
strong effect on the overall shape of the final nebula spectrum. 

The effect of this parameter on the gaseous emission can 
be seen, in particular, in the stronger Hi (eg. Balmer) edges as 
we move to lower ionization parameters. This is predominantly 
due to the lower electron temperatures which characterize the 
models with lower ionization parameter. The effect of dust can 
also be seen in several ways. For example, the greater overall 
extinction by dust, the greater flux in the UV arising from the 
input continuum scattered by dust, and the intensity of the dust 
far-IR emission bump. 

As we increase in ionization parameter, the fraction of UV 
flux absorbed by dust increases, and this combined with the 
increase in the ionizing intensity results in a greater total of 
IR emission by dust. It also results in hotter dust on average, 



which means the peak of the IR feature is shifted to shorter 
wavelengths and the flux in the mid-IR is greater. 

The effect of ionization parameter is more clearly seen in 
figure|3j which shows two models with the same input ionizing 
intensity, but different total pressures, giving different ioniz- 
ing parameters. The higher ionization parameter model (with 
Ptot/k - 10 7 K cirT 3 ) shows a stronger IR dust feature rel- 
ative to the lower ionization parameter model (P to t/k =* 10 s 
K cm -3 ). The dominance of dust absorption is most obvious 
however in the mid-IR to the optical-UV, where the log-scale 
emphasizes the separation in the models. This arises predomi- 
nantly due to the loss of ionizing photons to dust in the higher S 
model which results in a weaker free-free and free-bound con- 
tinuum. This increased absorption also produces more promi- 
nent FUV absorption. 

To demonstrate the relative features that arise due to pres- 
sure changes, figure |4] shows three model SEDs with the same 
initial ionization parameter (S = -0.66, log(f/o) = -1.883) but 
with different P tot . To account for the different input flux den- 
sities, the models have been scaled by 10.0, 1 .0, and 0. 1 for the 
models of P tot /k 10 6 , - 10 7 and 10 8 KcuT 3 respectively. 

Since they have the same ionization parameter, the gaseous 
line and continuum emission spectra of the models are very 
similar. However, the progression of the dust IR re-emission 
bump feature to shorter wavelengths as we increase in j tot is 
very marked. This progression can be explained by the increas- 
ing flux density experienced by the NLR cloud as the pressure 
increases to maintain the constant ionization parameter. The 
increasing incident flux density leads to higher average grain 
temperatures in the absorption region and thus stronger short 
wavelength emission. 

The UV spectra show that the same optical depth charac- 
terizes all three models. With the same fraction of the absorbed 
UV light heating the dust, the relative flux emitted by the dust 
in all three models is the same, as can be determined from fig- 
ure m 

When the models are allowed to continue to greater opti- 
cal depths, the lower flux densities mean that the dust emits 
at longer wavelengths. However, once the UV flux from the 
source is lower than the IR peak on a vF v plot, we don't get 
much more far-IR emission, rather the dust peak merely be- 
comes a little wider (see next section). Thus the dust is hotter 
at all optical depths for a greater incident flux density. 

Of special note is the presence of silicate emission features 
in the models with the highest incident flux densities. This in- 
dicates that in the most luminous AGN or most nuclear narrow 
line regions, such emission is possible. 

Figure|5]zooms in on the mid-IR region of the P to t/k 10 7 
K cirT 3 SEDs of figure|2] This figure reveals in more detail the 
silicate features, especially that at 10 fim. The 18 fim feature 
is less visible due to the shape of the IR continuum and its 
relative weakness. The figure also demonstrates more clearly 
the increasing maximum temperature achieved with higher in- 
cident fluxes, with the decreasing wavelength or the IR feature 
edge. 
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Fig. 2. Spectral energy distributions (SEDs) for NLR models with P lot /k =s 10 6 K cm -3 . The models decrease in ionization 
parameter from top to bottom, with log(£7o) = -0.4, —1.4, -2.2, -3.0 respectively. The SEDs presented here are purely nebular 
emission from the NLR. To obtain the full SED of the AGN, we would need to add these SEDs to the attenuated source SED, 
including any contributions from a Broad-Line Region (BLR) and/or accretion disk. 



5.1. The Effect of Column Depth 

One final parameter to consider in the determination of NLR 
dust emission is the total column density of the narrow line 
region. The previous models were all calculated to a total col- 
umn of N(Hi) = 10 2L5 crrT 2 , such that they were predomi- 
nantly ionization bounded. The column density affects the final 
line emission spectrum from the models as we either encom- 
pass (with increasing column density) or truncate (with lower 
column density) the various emission regions for each species. 

As far as the dust is concerned, the greater the column den- 
sity, the greater the optical depth of dust. At large column den- 
sities, most of the UV radiation field has been absorbed, and 
therefore the dust is at lower temperatures and emits weakly at 
longer wavelengths. Figure|6]shows the effect of increasing the 
final column depth in a model with pressure of P (ot /k 10 
Kern -3 and ionization parameter of H = -0.66. The column 
density is increased from 10 2ao cm -2 to 10 23 cm -2 in steps of 
0.5 dex. 

As expected, the low H i column depth models show less 
extinction by dust and have a narrower and hotter dust emis- 
sion bump. The highest column depth models have most of 
their optical absorbed by dust and show a much broader and 



stronger dust emission feature. The emission from the coolest 
dust arises predominantly from the largest optical depths in the 
models. This emission comes from the absorption of the low 
flux density visible and near-IR photons. 

The spectrum in the 10-30 fim is almost invariant with in- 
creasing optical depth until absorption becomes important, es- 
pecially in the strong 10 /imSilicate feature. This is because the 
flux here is entirely contributed by the hottest dust in the ion- 
ized gas or just behind the ionization front. Strictly speaking, 
as these models are 1-D the results become unrealistic for high 
column densities. More realistic 3-D treatments of the dust ra- 
diative transfer is required in these cases, especially as we are 
most likely viewing the NLR at various angles and do not in 
general have a direct line of sight through the obscuring mate- 
rial to the emitting edge of the NLR clouds. 

6. Discussion 

6.1. The Shape of the NLR IR spectrum 

It is only in recent times that the resolution of infrared tele- 
scopes has been high enough to begin to spatially resolve the 
emission from narrow line region, and only in the nearest of 
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Fig. 3. A comparison of two models having the same input ionizing flux density but differing total pressures of P tol /k st 10 7 K 
cm -3 (red curve) and P lot /k 10 s KcirT 3 (blue curve) demonstrating the differences due to ionization parameter and pressure. 



active galaxies JBock et allhoOOHGalliano et alll2005l) . In all 
other active galaxies it is difficult to distinguish the relative 
contributions of starbursts and of AGN in the nuclear regions. 
Nonetheless, it is important to try to do this, in order to under- 
stand the energetics of the nucleus as a whole, as well as the 
proposed AGN - starburst connection. 

A comparison of the model spectra shown in the previous 
section with th e recent high spatial resolution spectroscopy of 
NGC 1068 by Mas on et.alJ {2006) reveal several similarities. 
The observed NGC 1068 NLR spectra JMason et all [2006. Fig. 
2,), show that the wavelength of the spectral break between the 
nebula and dust continuum decreases with closer to the nu- 
cleus. This indicates hotter dust temperatures with increasing 
incident flux, exactly as in the models. The more nuclear spec- 
tra also show increasing silicate absorption due to increasing 
optical depth. This probably precludes the detection of any sil- 
icate emission from the inner NLR or nucleus itself. 

However the observed NLR spectra also reveal much flat- 
ter Mid-IR spectral slopes than obtained in most of the NLR 
models. The observations also indicate that the dust is hotter 
than found in our simple models discussed here. This suggests 
that there is some component of dust that the models are miss- 
ing. T his is supported by the measurements of Gratadour et al. 
(2006) who find constant high temperatures of ~500 K out to 



70pc from the nucleus of NGC 1068, where the NLR lies. This 
is much hotter th an can be predicted from simple NLR models 
jBarvainislfl987l) . 

The NLR of AGN most likely obtains a greater range of 
density and ionization than seen in our single cloud models 
which would soften this slope with the combination of both 
high and low temperature dust. More importantly it would also 
contain a mix of density-bounded (low Hn column density) 
and ionization bounded (low final H n fraction, as in the models 
presented here) clouds. This mix would allow the contribution 
of hot dust emission without significant absorption, as well as 
the stronger cool dust emission. This can also be accomplished 
though the combination of low covering factor hot dust, and 
large covering factor cool dust in a 3-D geometry. 

Smaller dust grains may provide another possible way for 
the models to reach higher temperatures. Currently, the average 
dust temperature of the front of the ionized cloud in the higher 
ionization P/k - 10 7 is ~250 K, and cooler temperatures as we 
step into the cloud. The smallest grains in the model, while hav- 
ing an average temperature lower than this, occasionally reach 
much higher temperatures due to stochastic effects, and pro- 
vide the hottest IR. They also tend to dominate the emission 
due to the grain distribution and larger surface area to mass. 
Thus, extending the distribution to even smaller sizes increases 
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Fig. 4. A comparison of three model SEDs with the same ionization parameter, log(t/o) = -1.883. The models have been scaled 
to remove the flux difference due to the different input flux densities (see text). As pressure decreases, the peak of the IR dust 
emission moves to longer wavelengths. 



the hottest dust emission, although it does not assist greatly in 
flattening the Mid-IR slope. 

Another possible explanation for the difference between the 
model and observed slope is alternative heating mechanisms 
not considered within these models, such as shocks. Although 
collisional heating is taken account of in the models, shock 
heating and destruction of dust grains can result in a much hot - 
ter IR spectrum, as considered in detail bv lContini et alJj2004l) . 
Even though emission line ratios indicate that photoionization 
is the dominant ionization mechanism of the NLR in most 
AGN, the correlation of jets with the NLR and high gas temper- 
atures and velocities observed within t he NLR a ll indicat e the 
shocks play a part in this region (see eg Almude na Prieto et all 
2005)). The combination of shock and radiation heating of dust 
may explain the observed slope and dust temperatures observed 
in the NLR of NGC 1068 and other AGN. 



6.2. The IR flux from the NLR 

While the contribution of the nucleus to the IR emission will 
presumably scale with the luminosity of the nucleus itself, the 
contribution of the NLR to the total IR emission of a galaxy 
will depend upon both the AGN luminosity and the NLR cov- 



ering factor. The models presented here allow, in principle, an 
estimation of the contribution of the NLR to the far-IR contin- 
uum. This exploits the coupling that exists between far-IR dust 
emission of the NLR and the emission lines generated within 
the NLR. 

Figure0shows the relationship between the strongest emis- 
sion line indicator for NLRs and AGN, the [O in] 500.7 nm line, 
and the IRAS 60 /urn flux (broadly equivalent to the the Spitzer 
MIPS 70 fim band). Note that the models have been scaled to 
remove the effect of source intensity, as an increase in source 
luminosity will result in both higher [O m] and 60 fim luminos- 
ity, when other parameters are kept constant. The ratio of the 
two fluxes shows some variation due to pressure, but more due 
to changing ionization parameter. At low ionization parameter, 
I[Oiir\ ~ 0.2 x I^Qfim, but at the highest ionization parameter 
I[onr\ ~ 0.04 x ho^m- However both pressure and ionization 

parameter can be estimated through the measure ment of other 

I 1 

emission lines such as [Sn] and Ha (see e.g. Groves etal., 
l2004bl) . reducing the possible range of the [Om]/60//m ratio. 

A potential problem in comparing the two fluxes is that the 
[Om] line may be heavily extinguished by dust along some 
lines of sight. In these cases, the theoretical ratio of the 60 yum 
flux to the [O m] 500.7 nm line can be used to provide an upper 
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Fig. 5. The Mid-IR spectra of the P lot /k ss 10 7 K cm 3 NLR models from[2] The models are shown here as log(vF y ) and are 
scaled to 1 at 5 //m. The 10 fim feature can be seen as a slight curve in the higher ionization models (top curves). 
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Fig. 7. The variation of the [Oin] 500.7 nm flux against the 
IRAS 60 fim flux for the NLR models. Shown are the P tot /k 
10 6 (afas/; c/of CMrve), ^ 10 7 (solid curve), and =; 10 8 (dashed 
curve), with the incident ionizing flux increasing from bottom 
to top along each curve as labelled. To remove the effect of 
source intensity variation the models have been scaled by 10.0, 
1.0, 0.1 for the P m /k * 10 6 , 10 7 and a 10 8 Kcm" 3 models 
respectively. 



estimate of the NLR luminosity. However, a better measure of 
the AGN luminosity absorbed by the NLR may be obtained 
using IR emission lines such as [Nev] 14.3 fim or [Oiv] 25.8 
yum. 

The 60 fim flux is also affected by the extent of the absorb- 
ing column of the NLR, with greater columns giving greater 
IR flux. A better measure would be to use the 25 //m band 
as it remains much less affected by the column depth, at least 
up to the point where Silicate absorption becomes important 
(~ 7V(Hi) > 10 2 , figure|6j. At these columns the IR emission 
lines will likewise be affected by dust extinction. The 25/im 
band is also where the NLR flux is most likely to dominate, ly- 
ing between the 12/im band, dominated by hot torus-like dust, 
and the FIR bands (60^m and 100/im bands) where star forma- 
tion is most likely to dominate. 

The [Nev] 14.3 fim flux is plotted against the IRAS 25 /mi 
flux in figure [8] The 25 fim continuum shows a much greater 
sensitivity to pressure due to the changing temperature of the 
hottest dust. The [Ne v] line is much more sensitive to the ion- 
ization parameter than the [O in] line. Together, these effects are 
expected to produce a scatter of ten in the correlation between 
the [Ne v] 14.3 fim and the IRAS 25 fjm fluxes. 



10 



B. Groves, M. Dopita, & R. Sutherland: The IR Emission from the NLR 



1 ,QPQP r 




0,1 



tO.Q 

wavelength Qwm) 



100.0 



Fig. 6. Spectral energy distributions (SEDs) for NLR models with P to t/k ==* 10 7 K cm 3 , S = -0.66, and varying in total Hi 
column. The models increase in Hicolumn depth from top to bottom (in the optical) from 10 200 cm -2 to 10 23() cm -2 in steps of 
0.5 dex. Note that the ionizing source is included to demonstrate the total absorption. 
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Fig. 8. The variation of the [Nev] 14.3 fim flux against the 
IRAS 25 fim flux for the NLR models. Curves and scaling are 
as in figure 



If we combine these lines we can get a good diagnostic for 
both the fraction of infrared emission arising from the NLR and 
the extinction experienced in the optical. Figure[9]demonstrates 



this with the ratio of [Ne v] 14.3 /mi/25 fim versus [O in] 500.7 
nm/25 fim. The first ratio gives an estimate of the NLR con- 
tribution to the 25 fim flux. If the observed value of this ratio 
is less than that given by the models, then it suggests that not 
all of the 25 fim flux arises from the NLR. However, the 25 
fim flux is strongly affected by the pressure of the models, or 
more precisely the incident flux density at a given ionization 
parameter. At a given ionization parameter, as the pressure and 
incident flux are increased, the [Ne v] 14.3/im line remains ap- 
proximately constant, at least for the lower densities. However 
with the increase in incident flux, the 25yum flux is increased, 
decreasing the [Nev]/25 fim ratio. Thus the pressure or inci- 
dent flux must be determined to give an accurate estimate of 
this ratio. In the IR this can be done using the density sensitive 
ratio of [Nev] 14.3/24.3 fim (see figure [Tot, thus connecting 
directly with this diagram. 

The [O m]/25fim ratio can also provide an estimate for the 
NLR contribution to the IR, but in addition gives an indica- 
tion of the obscuration of the NLR. As the NLR becomes more 
heavily obscured this ratio will be reduced by the loss of the 
[Oiii] line. 

This can be seen in the two observational points in this di- 
agram. NGC 1068 and NGC 4151 are two of the closest and 
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most studied active galaxies. Respectively these are quite of- 
ten used as the representative Type-2 and Type 1 AGN, corre- 
sponding to strongly obscured, edge-on, and lightly obscured, 
face-on AGN in the unified model. 

NGC 4151 appears to the lower left of the region occupied 
by the models. The observed mid-IR [Ne v] and [S n] lines give 
a density of ~ 10 3 cirT 3 (figure fTOt. indicating the ratios are 
less by 1 .0 to 1 .5 dex. This suggests that the NLR contributes 
only about 3 - 10% to the total 25jum emission, depending on 
the ionization state of the gas. 

NGC 1068 has a similar [Ne v]/25 jum ratio to NGC 4151, 
[Om]/25pm ratio of NGC 1068 appears less by ~ 1 dex com- 
pared to NGC 4151. While the higher density of NGC 1068 
indicated by the mid-IR lines (~ 10 4 cm 4 , figure ITUl and a 
differences in the average ionization parameter can possible ex- 
plain this difference, some part of the difference may also arise 
from the attenuation of the [O in] flux. This would fit within 
the Type-1 and 2 paradigm and is supported by the NGC 1068 
emission line asymmetries (Lutz et al., 2000) and estimates of 
low reddening in NGC 4151 iKriss et allll995l) . 

Assuming a higher density of ~ 10 4 cirT 3 and higher ion- 
ization parameter figure|9]suggests that the NLR of NGC 1068 
has a similar contribution of about 3 - 10% to the total 25/mi 
emission. 
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Fig. 9. Diagnostic [Nev] 14.3 fim/lRAS 25 fim versus [Om] 
500.7 nm/IRAS 25 fim ratio diagram. Each curve shows a dif- 
ferent pressure, decreasing from left to right. The direction of 
increasing ionization parameter for the curves is indicated by 
the values of Iq marked next to the P/k = 10 6 model. The sym- 
bols show the observed ratios of NGC 1068 and NGC 4151 
(as marked), with measuremen t errors dominated by line flu x 
calibration of -20 %. Data from I Alexander et alJ dl999i EoO0): 
ISturmetalJ ll 2002l) . 



Obviously this is a low contribution, and means that there 
are other effects not considered within the simple models pre- 
sented here. 

One possible consideration is the effect of aperture. The 
[Ne v] emission line arises only from the AGN, and likewise 
the [O m] is dominated by the AGN contribution. By contrast, 
the IRAS aperture is large and can include the contribution 



of star formation in both NGC 1068 and NGC 4151. Recent 
high spatial resolution observations suggest that this contribu- 
tion is not large, appro ximately 10% or less jRock etallEoOol 
Radomski et al .,2003). This does not preclude any nuclear star 
formation contribution within the ~ 4" resolution of these ob- 
servations. 

Another possibility is that the NLR of these objects is a dif- 
ferent metallicity to that considered here. Most AG N NLR ap- 
pear to be best fit by metallicities of about ~ 2Z JGroves et all 
2004b). An increase in metallicity would increase the dust ab- 
sorption relative to the gas and hence the IR emission. The in- 
crease relative to the [Ne v] line is more uncertain as it is not 
depleted onto dust, but will be weaker due to the lower gas tem- 
perature in higher metallicity gas. In addition to metallicity dif- 
ferences, an increase of the dust to gas ratio through increased 
depletion could also increase the contribution of the NLR to the 
total IR, though it must be balanced with the line emission. 

However the most likely possibility is that there are other 
components of IR emission not contained within our simple 
models. The most obvious of these is the putative torus, or more 
precisely, the equatorial dusty gas that is needed to explain the 
unified model, and probably associated with the accreting gas. 
This gas must be of a high enough density that it has no as- 
sociated forbidden emission lines, such as [Om], are emitted 
and close enough to the nucleus that hot dust temperatures are 
achieved. 

In addition to this there may also be dusty coronal gas, 
within the NLR that is so hot that species such as O m have 
been ionized. While dust may have only a short lifetime within 
this region due to sputtering, the covering fraction of this gas 
may be large enough to contribute to the total IR emission. 

Connected with this is the possibility of shock dust heating. 
As mentioned in the previous section, shocks are known to ex- 
ist within AGN, and may contribute to the total IR emission. 
IContini et alJ J2004I) were able to explain the IR spectrum and 
flux of several AGN using shocks alone, suggesting that they 
may be a viable contribution to the total IR. 

One other source of dust heating may be nuclear star 
formation. Nuclear starformation is known to exist in sev- 
eral AGN through the presence of PAH emission features 
frmanishi & WadaU2004 e g), and will contribute to the total IR 
emission without contributing to the [Ne v] emission. Nuclear 
star formation may contribute to other emission lines however 
such as [Om], though the optical-UV lines may be heavily at- 
tenuated. 



6.3. Mid-IR Emission Line Diagnostics 

Any nuclear star formation will also produce emission lines, al- 
though not of as high ionization potential as Ne v. A potential 
means to separate the contribution to the IR originating from 
AGN and that coming from circum-nuclear star formation is 
therefore to examine IR line flux ratios. Diagnostic ratios from 
ISO, such as that of polycyclic aromatic hydrocarbon (PAH) 
feat ures against the hot dust con t inuum or strong emissio n lines 
(Lutz et all fl99i iGenzel etafl Il998t ILutz et all Eooll) . have 
also been considered for this purpose. Here, we restrict our ex- 
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amination to the IR emission lines found in both starforming 
galaxies and AGN. 

The choice of IR emission line diagnostics is similar to that 
in the optical regime. The chosen emission line ratios must 
be strong, and generally obtainable within the same spectral 
range. One important difference is that extinction and redden- 
ing of the lines is no longer a strong criterion for line choice. 
Se veral good diagnostic ra tios have previously be en suggested 
bv lAlexanderetalJ dl999l) and lSturm et ail J2002I) . 

Two such ratios discussed in previous works are the den- 
sity sensitive line ratios [Nev] 14.3/[Nev] 2 4.3 jum, and [S n 
18.7/[S n] 33.5 yum, as shown in figure[T()l(see Alexan der et al. 
1999, for a description of density sensitivity). These ratios can 
help remove some of the uncertainty seen in the previous dia- 
grams. In this figu re (and figure sll0landfT2l the observational 
data is from Stur m et alJ J2002), with relevant discussions on 
the precision of the data found within. 

Both ratios consist of strong IR lines and remove direct ion- 
ization and metallicity effects. As [Ne v] is such a high ioniza- 
tion species, it arises purely from the NLR, and in particular 
gives an indication of the densities in the hot, highly ionized 
regions. The [S n] ratio gives an idea of the densities in the 
cooler parts of the NLR, but can also be contaminated by emis- 
sion from nuclear star formation. 

Such a contribution may explain the offset between the ob- 
served AGN and the models, where the [Ne v] ratio suggests a 
higher density than the [S n] ratio for several objects. However, 
no coincidence with associated starbursts is seen (starred ob- 
jects). Another possibility is that the region of the NLR where 
the [S n] lines arise is at a lower density than the [Ne v] region. 
Like previous work, we use the [Ne v] ratio for our estimated 
density. 
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Fig. 10. Density sensitive [Nev] 14.3 yum/[Nev] 24.3 //m ver- 
sus [Sn] 18.7 jum/[Sn] 33.5 yum ratio diagram. Each curve 
shows a different pressure, increasing from left to right, as la- 
belled. The symbols show the observed line rati os of nearby 
AGN (as labelled) from the ISO-SWS survey of ISturm et all 
Triangles represent Seyfert 2 galaxies, while squares 
are Seyfert Is, with the star filled objects representing those 
with associated star formation within the aperture. The average 
ratio error of the AGN is -28%. 



To separate the nuclear starformation from AGN activity, 
the line ratio diagnostics generally depend upon lines that have 
either an ionization beyond that found in star forming regions, 
such as [Ne v] 14.3yum, or are sensitive to the extended partially 
ionized region found around AGN, like [S i 11] 34 yum. Such ex- 
amples are shown in ISturm et"al The problem with 
the high ionization potential lines is that they are very sensi- 
tive to the ionization conditions of the nebulae. Similarly, the 
lines arising in the partially ionized region are generally mat- 
ter bounded, and are sensitive to the total column depth of the 
emitting region. 

Figure^2 snows a diagram which partly avoids these prob- 
lems. It uses three close proximity, strong neon emission lines; 
[Nen] 12.8 ^m, [Nev] 14.3 fim and [Nem] 15.5 yizm. The 
[Nev]/[Nen] ratio is a strong indicator of the AGN contri- 
bution, as a nuclear starburst (SB) will only contribute to the 
[Ne n] line and thus only weakens this ratio. Typical NLR val- 
ues for this ratio are 3-10. However as seen from the model 
curves in figure ^2 this rat i° is verv sensitive to the ioniza- 
tion conditions, thus is degenerate between AGN-SB mixing 
and ionization. The second neon ratio, [Nem]/[Nen], helps to 
remove this degeneracy. As shown by the model curves, it is 
much less sensitive to the ionization state of the gas. While 
starbursts produce both these emission lines, the ratio is much 
lower than found in AGN (depending upon the metallicity of 
the SB). 

The usefulness of this diagram can be seen in the distri- 
bution of the observed AGN, with all objects with a known 
starburst component lying to the lower left of the diagram. A 
possible empirical dividing line for starburst affected AGN then 
suggested by this diagram is 



[NeIII]/[NeII] < 0.9([NeV]/[NeII]) - 0.35. 



(5) 



A more theoretical mixing diagram will be presented in fu- 
ture work. 

One possible source of error in this diagram in stronger 
starbursting galaxies, is the proximity of the neon emission 
lines to the PAH emission bands, which may make mea sure- 
ment of these lines difficult (see e.g. Weedm an et all l2005). 

6.4. Distinguishing the NLR and Dusty Torus in IR 
Emission 

While the line diagnostic ratio diagrams are a good way to sep- 
arate the contribution of starformation and AGN to a spectrum, 
separating the contribution of the NLR and the putative inner, 
dusty torus is more difficult. The only possible way is to use 
line-continuum ratios, either as equivalent widths or line-band 
ratios. Figure|9]shows one possible diagnostic diagram, though 
this requires the subtraction of any star formation contribution 
to the 25 /jm flux to properly determine the torus contribution. 

A possible alternative to this is shown in figure[2] It uses a 
strong, high ionization emission line, in this case [Oiv] 25.9 
yum, against the 25 yum and 12 yum IRAS fluxes. This dia- 
gram is similar to the IRAS colour-colour diagrams (see eg 
iDopita et all fl998). but includes the Mid-IR emission line to 
account for the IR emission from the NLR. We expect the dusty 
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Fig. 11. Ionization source diagnostic diagram of [Nev] 14.3 
//m/[Nen] 12.8 //m versus [Nem] 15.5 /mi/[Neii] 12.8 /mi. 
The direction of increasing ionization parameter for the model 
curves is indicated by the values of Iq at each model point. The 
observed AGN ratios are marked as in figure[l0| with the same 
average measurement errors of 28%. Any contribution of star- 
formation will move objects to the lower left of the diagram. 



torus to be significantly hotter than the NLR, and thus show 
stronger emission at shorter wavelengths. Therefore the emis- 
sion at 12 /mi should be less correlated with the [O iv] emission 
line than the 25 /mi flux. 

The observations support this, with the observed AGN ap- 
pearing 0.5 to 1 dex lower than the models in [O iv]/25 /mi, but 
approximately 1 to 2 dex lower in the [O iv]/25 //m ratio. This 
diagram also gives an indication of a much greater contribution 
of the NLR to the IR flux, and reveals the need to use several 
line ratios to get a good indication of the NLR contribution. 
Depending on the density assumed, the [O iv]/25 /mi ratio in- 
dicates a NLR contribution of 10-50% to the 25/rni flux, and 
a weaker contribution to the 12/mi flux. Determining the exact 
contribution of the NLR to the total IR will depend upon both 
the shape of the torus IR feature and the diffuse IR emission. 

A NLR contribution of 10-20% to the mid-IR flux is not 
unreasonable and fits in with observations. High resolution ob- 
servations also show the dominance of the torus, with sub- 
arcsecond images of NGC 1068 indicating that around 75% 
of the Mid-IR emis sion arises fr om the unresolved core where 
the torus lies feocket alll2000l) . 

In terms of energetics, the contribution of the NLR de- 
pends critically to the covering fraction of the NLR which is 
an unknown quantity. If dust is ignored, a H/3 line emission 
analysis o f nearby Seyferts give s a covering fraction of only 
1% -4% JNetzer & Laorl Il993l) . However, when dust is in- 
clud ed it com petes with Hydrogen for the ionizing photons 
(Do pita et al.ll2002l) . and thus much greater covering fractions 
are possible. Likewise, observations of nearby Seyfert galaxies 
appear to suppo rt a cov ering fraction greater than a few percent 
JNetzer & Laorlll993l) . 

Another indicator that the NLR may have a large covering 
fraction is the observations of UV or "associated" absorbers. 
These absorption systems are believed to have a covering frac- 
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Fig. 12. IR line-continuum diagnostic diagram of [Oiv] 25.9 
/mi/ IRAS 12/mi band versus [O iv] 25.9 /mi/ IRAS 25/mi band. 
Observed AGN ratios are marked as in figure^)! The IR contri- 
bution from a dusty torus will move objects lower and strongly 
to the le ft. Averagee measu rement errors are approximately 
20% (see lSturmet alll2002i) . 



tion > 50% JCrenshaw et all 120031) . and several observation 
indicate a possible connection between these and the NLR 
JCecil et all l2002t ICrenshaw & Kraemerl 120051) . While these 
systems do not have the column density observed in the NLR, 
figure [6] shows that such systems can still contribute to the 
IR emission if they contain dust. Thus a NLR contribution of 
~ 10% to the total IR flux in AGN is not unfeasible. 

As a last point, while these models do not represent in 
any form the putative torus around the AGN from the unified 
model, this work can be extended to this region, albeit with- 
out the consideration of geometry. As we increase the density 
and incident flux onto the models we start to reach the possi- 
ble physical conditions within the torus, as the appearance of 
silicate emission in figure[5]shows. 

7. Conclusion 

We have explored within this work a series of Dusty Narrow 
Line Region (NLR) models, covering a range of realistic pres- 
sures and ionization parameters, with the aim of investigating 
the Infrared emission from dust in these regions. 

The presented spectral energy distributions (SEDs) for 
these models reveal the possible range of the IR emission from 
the NLR, and the correlation of this emission with the ioniza- 
tion parameter and the intensity of the ionizing radiation. An 
increase in either of these parameters will increase both the to- 
tal IR flux from the dust and the peak temperature of the dust. 

Comparisons with observed NLR IR spectra show that sim- 
ple 1-D photoionization models are not enough to represent the 
complexity of the NLR. Combinations of different NLR clouds 
and NLR dynamics may be needed to represent a more accurate 
model of the NLR IR emission. 

While examining the spectral energy distributions of these 
models gives an idea of both the overall emission from the NLR 
of active galaxies, the greatest benefit of these models is that 
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they give a direct connection between the emission by dust 
and the line emission from the photoionized gas. A selection 
of correlations between line and dust emission are presented 
here to give possible diagnostics for the IR emission from AGN 
and galaxies. In addition we have explored further possible IR 
emission line diagnostics to assist in the determination of the 
NLR physical conditions, and to help distinguish the contribu- 
tion that starformation makes to the total spectrum. 

Comparisons of the models with the nearby AGN, NGC 
1068 and NGC 4151, appear to show that only -10% of the 
25/mi emission arises from the [O m] line emitting region. This 
means a significant fraction must arise from either the putative 
torus, nuclear star formation, coronal-shock excited region, or 
some combination thereof. 

While the fraction the NLR contributes to the overall in- 
frared continuum of an active galaxy may be small, it is impor- 
tant to consider this component when trying to understand the 
full spectral energy distribution of AGN. The final aim will be 
to connect the emission from the NLR with what is occurring 
in the central engine and thus connect the emission lines with 
the total IR output of the AGN. 
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